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Abstract 

The 4/ — 5d transition rates for rare-earth ions in crystals can be calculated with an effec- 
tive transition operator acting between model 4/^ and 4f^~^5d states calculated with effective 
Hamiltonian, such as semi-empirical crystal Hamiltonian. The difference of the effective transition 
operator from the original transition operator is the corrections due to mixing in transition initial 
and final states of excited configurations from both the center ion and the ligand ions. These 
corrections are calculated using many-body perturbation theory. For free ions, there are impor- 
tant one-body and two-body corrections. The one-body correction is proportional to the original 
electric dipole operator with magnitude of approximately 40% of the uncorrected electric dipole 
moment. Its effect is equivalent to scaling down the radial integral {5d\ r |4/), to about 60% of the 
uncorrected HF value. The two-body correction has magnitude of approximately 25% relative to 
the uncorrected electric dipole moment. For ions in crystals, there is an additional one-body cor- 
rection due to ligand polarization, whose magnitude is shown to be about 10% of the uncorrected 
electric dipole moment. 
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I. INTRODUCTION 



In the recent years, optical spectroscopy of lanthanides ions in crystals involving 4/ — 5d 
transitions has been widely studied. These transitions are electric dipole parity allowed, and 
therefore suitable for efficient absorption of VUV radiation from the noble gas discharge in 
mercury-free lamps and plasma display panels-i^^^ A well-known application of the strong 
4/ — 5d absorption of lanthanides is in blue lamp phosphors BaMgAlioOiy: Eu^^. Potential 
applications for which VUV 5d — 4/ emission can be used are tunable VUV lasers and 
scintillator materials. 

The transition rates are due to electric dipole moment contribution. The electric dipole 
matrix elements between 4/ and 5d orbitals have been generally expected to dominate the 
electric dipole moment, and most of the calculationa^*^^ we are aware of consider only 
this zeroth-order contribution, which we shall refer to as the uncorrected electric dipole mo- 
ment in the paper. Recently, various detailed comparisons between theory and experiment 
for free ions have shown that the zeroth-order calculations considerably over-estimated the 
5d 4/ spontaneous emission rates and configurational mixing needs to be included.— i^iiS 
Using the pseudorelativistic Hartree-Fock (HFR) approach of Cowan^^, a configuration in- 
teraction (CI) calculation for free-ion Ce^+ includes both 12 odd-parity and 12 even-parity 
configurations gave radiative lifetimes agrees to experiment,— while a single configuration 
calculation gives uncorrected radial integral {5d\ r |4/) = 0.0435nm, which underestimates 
the lifetime, i. e., overestimates the rates, by a factor of about 3. For rare-earth ions in 
crystals, the calculation of electric-dipole allowed 4/ — 4/ two-photon absorption rates has 
shown that some contributions due to ligand excitations are very important.— Since the 
ligand polarization contributions to two-photon transition moment are mainly due to mod- 
ification of one-photon transition moment between 4/^ and Af^'^Ad states, they are also 
potentially important to 4/ — 5d one-photon transitions. 

Direct CI calculations similar to the one for Ce^"*" are less feasible for rare-earth ions 
with near half-filled 4/ shell due to the large dimension of state spaces and become even 
more impractical when ligands orbitals are included, as in the cases of quantum chemical 
cluster calculations. A systematical method to take all those contributions into account in a 
reduced space is effective operator method . ^^i^^i^^ In this method, the full Hilbert space time- 
independent Hamiltonian H is transformed into an effective Hamiltonian ifefr, which acts on 
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the reduced space (referred to as model space) and gives upon diagonalization a set of exact 
eigenvalues and model space eigenvectors. For a time- independent operator O, an effective 
operator Ogfr may be introduced that gives the same matrix elements between the model 
space eigenvectors of H^s as those of the original operator O between the corresponding true 
eigenvectors of H. Effective Hamiltonians H^g and transition operators O^s are often con- 
structed by many-body perturbation theory (MBPT) with order by order approximationi^ 
and then represented by connected diagrams similar to Feynman diagrams. With those dia- 
grams, perturbation calculations involving great number of many-body intermediate states 
transform into calculations involving summations over one- and two-particle states. Hence 
the calculations are usually greatly simplified and the results are much easy to interpret 
since we can find out transition mechanism from diagrams with important contributions. 
The calculations can also be easily used together will semi-empirical crystal Hamiltonian 
for 4/^ configuration^^ and its adaption to 4:f^~^5d configuration,— and with more general 
ligand field approaches. 

In this paper we do a general many-body perturbation calculation for the 4/ — 5d effective 
transition operator. The zeroth and first-order contributions to the effective electric dipole 
operator are presented in Sec. II in terms of diagrams and algebra expressions. In Sec. 
Ill we make further approximations to select out important contributions and discuss the 
possible consequences to 4/ — 5d one-photon and 4/ 4/ two-photon transitions. 



II. GENERAL MANY-BODY PERTURBATIVE EXPANSION FOR EFFECTIVE 
OPERATOR AND DIAGRAMMATIC REPRESENTATION 

We partition the Hamiltonian H for the center-ion-ligand system as usual into three 
terms, 

H = H0 + V1 + V2 (1) 

where Hq is a zeroth order model Hamiltonian, Vi and V2 are the part of one-body spin-orbit 
and two-body Coulomb interactions, respectively, not included in Hq. Usually one choose 
Hq to be spherical and spin independent, here we retain spin-independence but do not limit 
Ho to be spherical, so that if necessary the strong crystal-field interactions for 5d electron 
in crystals can be included in Hq, as in the calculation of 4/ — 4/ two-photon absorption 
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spectra in Eu^^:CaF2.— In general, the effect of crystal-field interactions to f — d transitions 
is mainly through the effect to eigenstates of the center ion. 

The eigenstates of Hq are antisymmetric product of one-particle states from the center 
ion and ligand ions. The model space is chosen to contain the product of both 4/^ and 
4f^~^5d states with states of completely filled shells from the ligands. 

The effective Hamiltonians can be calculated with a many-body perturbation expansion.— 
However, for 4/^ energy levels of rare-earth ions, the phenomenological crystal-field 
Hamiltoniani^ with adjustable parameters turns to produce better energy levels. The phe- 
nomenological crystal-field Hamiltonian has also been adapted to the calculation of 4/^~^5(i 
energy levels.— Since the eigenvectors of the phenomenological Hamiltonian can be used as 
the model space eigenvectors of the Hermitian effective Hamiltonian, we need only construct 
the Hermitian effective transition operator. 

We consider transition due to electric dipole mechanism only. The transition operator is: 



In addition to the zeroth-order contribution included in most calculations, here we include 
also the first-order contribution. Using our latest results on effective operators and modified 
diagrams we get the zeroth- and first-order diagrams for the effective transition operator 
as shown in Fig^ 

Since the effective operator acts between eigenvectors of effective Hamiltonian, which are 
linear combinations of bases in model space, the incoming lines and outgoing lines in FigQ 
can only be 4/ orbitals and 5d orbitals. The internal lines can be any orbitals of the center 
ion and ligands, as long as at any horizontal level at least one of them is a core or unoccupied 
orbital. 

The evaluation rules have been summarized previously by Duan et al!^^ For these diagrams 
in Fig. ^ the contribution from each diagram is proportional to the product of matrix 
elements of all vertexes divided by the net outflow energy (outflow energy taking away 
inflow energy) of the V vertex in the diagram. The algebraic expressions are as follows: 
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where m, n, p, g are valence orbitals, i.e., f or d orbitals, i can be any core or unoccupied 
orbital of the center ion or ligand ions, u can be a valence or unoccupied orbital of the center 
ion or an unoccupied orbital of the ligand ions, a is a core orbital of the center ion or ligand 
ions. Note that the two-body matrix elements are antisymmetric, i.e. 



iV2)ab,cd = {ab\ V2 \cd) - {ah\ \dc) . 



(8) 



III. FURTHER APPROXIMATIONS AND IMPORTANT CONTRIBUTIONS 

We make the following two approximations: (1) we follow the usual approximation by 
neglecting integrals contain overlap of 4/ orbitals of the center ion with ligand orbitals; (2) we 
neglect the extremely weak spin-orbit interactions between two orbitals nl and n'l {n ^ n').— 
With these approximations, the one-body contributions and two-body contributions can be 
greatly simplified and the strength can be calculated using structure data and HFR results 
for free ions. 



A. One-body contributions 

The zeroth-order term is the electric dipole between 4/ and bd orbitals, denoted by 
Fig^a). The first-order diagrams Fig^b,c) can be neglected with the second approxima- 
tions above. This leaves the only one-body contributions from FigCJd). Under the first 
approximation above, there are two kind of contributions from Fig^d): (1) the contribu- 
tion with the two internal lines are orbitals of the center ion, which does not change much 
from free ion to ions in crystals; (2) the contribution with the two internal lines are orbitals 
of the ligands, which is usually denoted as "dynamical coupling" or "ligand polarization" 
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contribution. The contribution due to ligand-center ion overlap are neglected under the first 
approximation since it contains 4/-ligand overlap. 

The first contribution. For the case of free ions, since both the Coulomb interaction 
V2 and the energy denominator in the expression for Fig^d) are spherically symmetric 
and spin-independent, the resulting interaction for FigC^d) can only be one-body spin- 
independent rank-1 spherical tensor, which is proportional to the electric dipole between 4/ 
and 5d orbitals with opposite sign. This provides a solid ground for calculations adopting 
an effective radial matrix element {5d\ r |4/)gjj scaled down from {5d\ r |4/). This contribu- 
tion is dominated by the terms with the pair of internal lines {5p, 5d), which have small 
energy denominators and large matrix elements. After complex manipulation of coupling 
and recoupling coefficients, we get the effective operator in terms of radial integrals and 
(re) coupling coefficients as follows: 
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where C. C. denotes the complex conjugative term, {ii, 12) is a pair of single-electron radial 
wavefunction indices, with one being core orbitals and the other valence or unoccupied 
orbitals. The leading terms are 
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(14) 

35(5rf|r|4/)(2e5rf-e5p-e4/) ^ ' 

The values of 5 calculated from the leading terms are 0.3571 and 0.3502 for Ce^+ and Pr^+ 

respectively. Note that 5 or Defr does not depend on the crystal environment. 

For free ion Ce^^, since there is only one active electron in model space 4/-|-5(i, the whole 

effective operator is strictly a one-body operator, which is approximately proportionally 
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to the original electric dipole D. The lifetime values from recent measurement and CI 
calculatioi>i2, can be used to work out the effective {5d\ r |4/)gjj ^ 0.025nm. Compared to 
the uncorrected HFR value {5d\ r 14/)^^ = 0.0435nm,^^ this gives an experimental 6 value 
for Ce3+, which is (1 - {5d\ r |4/)^g / {5d\ r |4/)) ^ 0.43. 

The second contribution. Neglecting the overlap between orbitals of the center ion 
and those of the ligand ions, we can do a bipolar expansion to the coulomb interaction. We 
obtain the same result as has been given earlier by Reid and Richardsoii^^ in calculation of 
two-photon transitions: 



Deff((rf),2) = (-e) Y: {5d\r'-'\Af)[ik)i2k-l)i2k + l)/3]'/'x 

fc=2,4,6 

E«^(e5d - e,f)Rl^'^'\c\L)C'-'^^^]^ (15) 

L 

= (-e(5rf|r|4/)) V2kTl{A\at,d,f)'-')^ + C. C, (16) 

fc=2,4,6 



_ /t M Ak - 1) |4/) (5rf||C^-^||4/) a,{e,, - e,f)C'^{L) 
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where L labels ligands, i labels the valence rare-earth electrons, and 



«lH = ^ E E m r, p ( + i ) (18) 

is the isotropic polarizability of ligand L. In the summation, q is over the component of 



and c and u are over core and unoccupied orbitals of ligand L, respectively. Here we 
introduce dimensionless coefficients { k = 2,4,6, q = —k, —k + 1, ■ ■ ■ , k) with associated 
operators properly normalized so that their values reflect the relative magnitudes of the 
contributions compared to uncorrected electric dipole moment. 

For Ce^^: CaF2, we use the structure data for CaF2, setting ~ 10~^nm^ for and 



Qg3+ ^Yee-ion data^ito obtain the A'^ values. The nonzero coefficients are Aq = y 14/5^4^4 = 
-0.08 and A^ = -^I2/1A\^ = 0.0064. The case for 4/ - 5d transitions of Gd^+iLaFg is 
estimated to be similar. Therefore, corrections due to ligand polarization is not so important 
for 4/ — 5d transitions as the corrections due to the excited states of the center ion. The 
reason that the ligand polarization contributions are important for ^3^/2 — > ^Ij 4f —>■ 4/ 
two-photon absorption is due to superposition of contributions from different intermediate 
states: the transitions are between states with main components satisfying AL = 6, most 
part of the contribution due to the 4/ — 5d rank-one uncorrected electric dipole contributions 



cancels, but the rank-3 and rank-5 ligand polarization contributions do not cancel and 
become important. 



B. Two-body contributions 

The two diagrams Fig^e) and Fig^f) are complex conjugation of each other. To give 
a nonzero contribution, the single internal line can only be an orbital of the center ion due 
to the negligible overlap between 4/ and ligand orbitals. These two diagrams are likely 
to be dominated by the two terms with 5p internal lines, which have energy denominators 
and matrix elements comparable to the important term in Fig^d). The two terms with 
5p internal lines, which are complex conjugation of each other, can be written as (with the 
complex conjugate term neglected) 
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where i^ = 2,4, Q = K — 1, K, and K + 1 and the dimensionless coefficients 
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The operators associated with d{K, Q) are appropriated normalized so that the relative 
strengths can be reflected by these d{K, Q) coefficients. The d{K, Q) values for Pr^+ are given 
in Table|l] These values show that the corrections due to two-body effective operators have a 
magnitude of 25% and can be important for rare-earth ions with two or more electrons in 4/ 
and 5c? shells. Since in general matrix elements of two-body operators are not proportional 
to those of one-body operators, the effect due to two-body effective operators cannot be 
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fully accounted for by scaling the electric dipole operator. 

All the different contributions to 4/ — 5d electric dipole transitions we considered and 
their relative magnitude are listed in Table IHl The uncorrected electric dipole term and 
one-body free-ion term are proportional to each other and form the main contribution to 
4/ — 5d transitions, which are equivalent to scaling the 4/ — 5d electric dipole radial integral 
with a factor 60% (= 1 — 40%). We call this corrected one-body free- ion term. The next 
important contribution is the two-body free-ion term, which has a relative magnitude of 
about 25% of the uncorrected electric dipole term and 40% (~ 25%/60%) of the corrected 
one-body free-ion term. This term has the same rank as corrected one-body free-ion term 
(rank-1), but is two-body and cannot be fully accounted for by scaling the 4/ — 5d electric 
dipole radial integral. The ligand-polarization term is about 10% or the uncorrected electric 
dipole and 15% of corrected one-body free-ion term. This term is less important than two- 
body free-ion term for 4/ — 5d transitions but the high rank (rank-3 and rank-5) part can 
be important for 4/ — 4/ two-photon transitions with AL > 2. 

IV. SUMMARY 

In summary, we have presented a first-order many-body perturbation calculation of the 
effective operator for 4/ — 5d transitions to account for the main corrections to the transition 
rates due to mixing of other configurations of the center ion and ligand orbitals in 4/^ 
and 4:f^~^5d states. Further approximations are made to select out possible important 
corrections to the usual zeroth-order 4/ — 5d electric dipole operator. First-order one-body 
contributions due to excited states of the center ion is shown equivalent to scale down the 
electric dipole radial integral {5d\ r |4/), which, for Ce^"^ and Pr^+, has almost reduced the 
effective value to about 60% of the uncorrected values. For rare-earth ions with two or more 
electrons in 4/ and 5d open shells, first-order two-body contributions also have important 
contributions, which cannot be accounted for by scaling the electric dipole radial integral 
{5d\ r |4/). The magnitudes have been calculated to be about 25% relative to the uncorrected 
electric dipole moment for Pr^+. Contributions due to ligand polarization are rederived, 
which are the same as given earlier by Reid and Richardson.— The magnitude is about 10% 
relative to the uncorrected electric dipole moment for Ce^+:CaF2. The contribution due to 
ligand-center ion overlap is negligible, since it contains the overlap between 4/ and ligand 
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orbitals. 
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Figures 



FIG. 1: Zeroth and First order diagrams for the effective electric dipole transition operator, where 
the free lines are 4/ and 5d orbitals, the internal lines in diagrams (b), (c), (e) and (d) are core 
or unoccupied orbitals, one of the two internal lines in (d) is core orbital and the other can be 



unoccupied or valence orbitals. The rules to evaluate these diagrams can be found in Ref 



2C 



12 



Tables 



TABLE I: The values of dimensionless coefficients for Pr^+ two-body transition effective operators, 
where c/tot = iJ2d{K, Q)^)^/^ denote the total magnitude of two-body interactions. 

TABLE IL List of different contributions to 4/ — 5d transitions and their relative magnitude 
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Figure d Duan and Reid, Journal of Chemical Physics 
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Table U Duan and Reid, Journal of Chemical Physics 



d{21) 


d(22) 


d(23) 


d{43) 


d{U) 


d{45) 




-0.147 


-0.155 


-0.113 


-0.015 


-0.035 


-0.063 


0.253 
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Table m Duan and Reid, Journal of Chemical Physics 



uncorrected 


one-body free- ion 


two-body free-ion 


ligand-polarization 


1 


(-)40% 


25% 


10% 
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